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EXPONENTIAL GROWTH RATES IN A TYPED 
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We study the high temperature phase of a family of typed branch- 
ing diffusions initially studied in [Asterisque 236 (1996) 133-154] and 
[Lecture Notes in Math. 1729 (2000) 239-256 Springer, Berlin]. The 
primary aim is to establish some almost-sure limit results for the long- 
term behavior of this particle system, namely the speed at which the 
population of particles colonizes both space and type dimensions, as 
well as the rate at which the population grows within this asymp- 
totic shape. Our approach will include identification of an explicit 
two-phase mechanism by which particles can build up in sufficient 
numbers with spatial positions near — ft and type positions near 
Ky/i at large times t. The proofs involve the application of a variety 
of martingale techniques — most importantly a "spine" construction 
involving a change of measure with an additive martingale. In ad- 
dition to the model's intrinsic interest, the methodologies presented 
contain ideas that will adapt to other branching settings. We also 
briefly discuss applications to traveling wave solutions of an associ- 
ated reaction-diffusion equation. 

1. Introduction. In this article we will consider a certain family of typed 
branching diffusions that have particles which move (independently of each 
other) in space according to a Brownian motion with variance controlled by 
the particle's type process. The type of each particle evolves as an Ornstein- 
Uhlenbeck process and this type also controls the rate at which births occur. 
The particular form of this model permits many explicit calculations, but 
throughout we will strive to develop techniques that rely on general prin- 
ciples as much as possible, so they might readily adapt to other situations. 
This model was previously considered in [12, 13]; these papers form essential 
foundations for this work, although we will recall various results as necessary. 
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We will make some significant applications of the spine theory for branch- 
ing processes. Inspired by the series of papers Lyons, Pemantle and Peres 
[19], Lyons [18] and Kurtz, Lyons, Pemantle and Peres [16], spine techniques 
have been instrumental in recent years in providing intuitive and elegant 
proofs of many important classical and new results in the theory of branch- 
ing processes. In this article we use the recent reformulation of the spine 
method presented in [8], which follows in similar spirit to the branching 
Brownian motion study of Kyprianou [17]. For a selection of other applica- 
tions of spine techniques, for example, see [1, 6, 7, 23] and references therein. 

1.1. The branching model. We define Nt to be the set of particles alive 
at time t > 0. For a particle u E Nt, X u (t) E K is its spatial position, and 
Y u (t) S 1 is the type of u. We will label offspring using the Ulam-Harris 
convention where, for example, if u = 021 then particle u is the first child 
of the second child of the initial ancestor, and we will write v > u if particle 
v is a descendant of particle u. The configuration of the branching diffusion 
at time t is given by the point process X t := {(X u (t), Y u (t)) :u E N t }. 

A particle's type evolves as an Ornstein-Uhlenbeck process with an invari- 
ant measure given by the standard normal density <j)(y) and an associated 
differential operator (generator) 



where 9 > is considered to be the temperature of the system. The spatial 
motion of a particle of type y is a driftless Brownian motion on ]R with 
variance 



Each offspring inherits its parent's current type and spatial position, and 
then moves off independently of all others. We use P x ' y and E x ' y with x, y E 
M to represent probability and expectation when the Markov process starts 
with a single particle at position (x,y). 

We will find the almost sure rate of exponential growth, D(p/,K), of par- 
ticles which are found simultaneously with spatial positions near — *yt and 
type positions near ny/i at large times t. From this we can deduce the speed 
of extremal particles and hence the asymptotic shape of the particle system. 
The main effort is required in identifying D(j, k) as the almost sure limit of 




A(y):=ay 2 , 



where a > 0. 



A particle of type y particle is replaced by two offspring at a rate 
R(y) := ry 2 + p, where r, p > 0. 



{X u (t)<-ft;Y u (t)>KV~t}- 



ueN t 
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In particular, the convergence properties of two different families of addi- 
tive martingales associated with the branching diffusion will lead directly to 
the spatial exponential growth rates and an upper bound on the space-type 
growth. For the remaining lower bound, we describe an explicit two-phase 
mechanism for amassing the required number of particles with prescribed 
space-type positions. The first phase involves building up an "excess" num- 
ber of particles, each covering a certain proportion of the required spatial 
distance. During their second phase, enough of these particles must succeed 
in making a difficult and rapid ascent into the required position. The latter 
phase is proved using an intuitive change of measure technique that induces 
a spine construction. 

The family of models we are considering is specific but nevertheless have 
some features of fundamental significance that motivate the choices for Qg, 
R and A. If the spatial motion is ignored, we have investigated a binary 
branching Ornstein-Uhlenbeck process in a quadratic breeding potential. In 
contrast, Enderle and Hering [5] considered a branching Ornstein-Uhlenbeck 
with constant branching rate but random offspring distribution. A quadratic 
breeding potential is a critical rate for explosions in the population of parti- 
cles. In a branching Brownian motion on K with binary splitting occurring 
at rate x p at position x, the population will explode almost surely in finite 
time if p > 2, whereas for p = 2 the expected number of particles explodes 
while the total population remains finite for all time with probability 1 (see 
[15], Chapter 5.12). The Ornstein-Uhlenbeck process is not only a canoni- 
cal ergodic diffusion, but this type-motion has exactly the right drift to help 
counteract the quadratic breeding rate. For high temperatures, 8 > 8r, there 
is a sufficiently strong mean-reversion in the type processes to ensure that 
the expected total population size does not blow up; but for temperatures 
9 < 8r, the quadratic breeding overpowers the pull toward the origin, the 
expected population blows up in a finite time and particles behave very dif- 
ferently. Throughout this paper we consider only high temperatures 6 > 8r, 
deferring the low and critical temperature regimes to future work. Given 
other choices, the quadratic spatial diffusion coefficient now becomes very 
natural, enabling us to find explicit families of (fundamental) additive mar- 
tingales since the linearized traveling-wave equation can be linked to the 
classical harmonic oscillator equations from physics. The binary branching 
mechanism was taken for simplicity; in principle our approach could extend 
to general offspring distributions, although new features would arise from 
possible extinctions and necessary offspring moment conditions. All these 
choices make the models rich in structure, possessing some very challenging 
features whilst remaining sufficiently tractable. 
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1.2. Application to reaction-diffusion equations. Following in the foot- 
steps of McKean [20] , the solution of the reaction-diffusion equation 



with initial condition f(x,y) € [0, 1] for all x,y £ M, can be represented by 



Of great importance for reaction-diffusion equations are traveling-wave so- 
lutions (e.g., see [21]). In the present context, a solution to equation (1) of 
the form u(t,x,y) := w(x — ct,y) is said to be a traveling-wave of speed c, 
where w(x,y) solves the traveling-wave equation 



Fundamental to our study of the branching diffusion are two families of "ad- 
ditive" martingales, Z^(t) [defined at (6)], which are linked to the lineariza- 
tion of (1). When 6 > 8r, Harris and Williams [13] determined when is 
uniformly integrable (see Theorem 17) and then w\(x, y) := E x,y exp(— Z^~(oo)) 
yields a traveling wave of speed c^. This gives the existence of traveling 
waves for all speeds c greater than some threshold c{0) := inf c^. 

Furthermore, combining the McKean representation (2) with the almost- 
sure convergence result established in [12] (look ahead to Theorem 18) can 
give results on the attraction toward traveling waves from given initial data. 
For example, if — In/ (x, y) ~ e Xx g(y) uniformly in y as x — > oo for some suit- 
able g € L 2 ((p), the solution u(x,y) to (1) with initial conditions / satisfies 
u(t, x — c7t, y) — > w\(x + x, y) as t — ► oo, where x is some constant that can 
be determined from g. 

In future work we hope to develop the approach used for standard BBM 
and the FKPP equation in [11], and prove that traveling waves of a given 
speed c > c{6) are unique (up to translation) and that no traveling waves 
exist for speeds c < c{6). We anticipate that our new results on the growth 
rates of particles will aid in establishing some difficult estimates on the tail 
behavior of any traveling wave, and hence assist in proving the conjectured 
uniqueness. In addition, we expect our growth rate results will be essential 
in obtaining broader classes of initial conditions that are attracted toward 
traveling waves. In each of these problems, difficulties arise from the un- 
bounded type space where, for example, some control must be gained over 
the possible contributions to J2ueN t 1°§ f(X u (t) — ct,Y u (t)) from particles 
that have large type positions in addition to large spatial positions. 



(1) 





(3) 
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2. Main results. In this section, we will present our main results that 
identify the growth rates found within the branching diffusion. We will give 
an overview for our proofs, identifying the key ideas and techniques used, 
as well as introducing some intuition for the dominant behavior of particles 
that underpins our approach. 

2.1. Martingales. The principal tools used throughout this paper are 
two fundamental families of "additive" martingales, which were introduced 
in [13]. 

Before defining the martingales we give some key definitions. Let 



A, 



8r 



V 4a 

Let Asl, with the following convention which we always use for A: 

Amin < A < 0. 

Also, define 



(4) (i x :=ye(e-8r-4a\2), := \ ± 

(5) Ef:=p + 9^, cf:=-Ef/X. 

Will will occasionally write E^ as E^^X) in order to emphasize that E^ are 
really functions of A; the ± superscripts will always distinguish these from 
expectation operators. Note that A m j n is the point beyond which ji\ is no 
longer a real number. 

The martingales are and Z^, defined for A £ (A m i n ,0] as 

(6) Z±(t) := £ vi{Y u {t))e^- E t\ 

ueN t 

where v^(y) :=exp(^y 2 ) are strictly-positive eigenfunctions of the opera- 
tor 

Qe + \\ 2 A + R, 

with corresponding eigenvalues E^ < Et and A, R are the functions defined 
in Section 1.1. This operator is self-adjoint on L 2 ((j)) with the inner product 
(■,-)cj> where (f,g),p '■= J fg4>dy and <p is the standard normal density. Note 
that £ L 2 (</>), whereas ^ L 2 (c/)) so is not normalizable. 

The calculations of Section 3 make it easy to see these are martingales, 
and throughout the paper we will need a variety of martingale convergence 
results which are gathered together in Section 8. In particular, we will need 
to know precisely when Z^ is uniformly integrable with a strictly positive 
limit, some further strong convergence results for other closely related sums 
over particles (also identifying which particles contribute nontrivially to their 
limits), and the rate of convergence to zero of the Z^ martingales. 
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2.2. The asymptotic growth-rate of particles along spatial rays. As an 
essential initial step toward determining the growth rate of particles in the 
two-dimensional space-type domain, we first look at the growth rate of par- 
ticles in the spatial dimension only. 

For 7 > and C C R, define 

( 7 ) N t ( T ,C):= 1 {X^t)<- 1 t;Y u {t)dC}- 

ueN t 

The limit giving the expected rate of growth, 

lim t" 1 log£;(iV t (7;R)) 

t— >oo 

can be shown to exist and its value can be calculated to be 

A(7):= inf {E~ x + A 7 } 

(8) „ . 



= P+ l --J^(6-Sr){A^ + 6a). 

An outline for this expectation calculation is given in Section 3. 

It is now tempting to guess that the asymptotic speed of the spatially 
left-most particle, c(6), is given by 

c{6) := sup{ 7 : A( 7 ) > 0} 
(9) , 



2{2r + p) 2 

'2a[r + p + 



0-8r 

Recall that 5(9) = inf^ g (^ min ) cT/ is also the minimum threshold for traveling 
waves. In this particular situation, the guess that "expectation" and "almost 
sure" right-most particle speeds agree was first proved rigorously using a 
martingale change of measure technique in [13]. In this paper, we extend 
this connection and prove that the "expected" and "almost sure" rates of 
growth of particles with given speeds (Theorem 1) and given space-type 
locations (Theorem 3) agree. 

Theorem 1. Let 7 > and yo<yi- Under each P x ' y law, the limit 
D(j):= lim t- 1 log^(7;[?/o,yi]) 



exists almost surely and is given by 
D( 7 ) = 



A( 7 ), if0< 7 <c(9), 
-oo, if 7 >c(0). 
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Note that symmetry in the process means there is a corresponding re- 
sult for particles with spatial velocities greater than +7 (corresponding to 
positive A values). We may occasionally make use of such process symme- 
tries without further comment. Then, since N^jM) is integer valued, the 
asymptotic speed of the right-most particle follows immediately: 

Corollary 2. Almost surely, 

lim t^ 1 sup{X u (t) :u£N t } = c(0). 

t— >oo 

This spatial growth rate result is proved in Section 10 using the martin- 
gale results from Section 8. In fact, it is very easy to obtain the upper bound 
by first dominating the indicator function with exponentials to reveal that 
N t (j;R) < exp{(E7 + Xj)t}Z7 (t), recalling that Z7 is a convergent mar- 
tingale, and then optimizing over the choice of A. For the lower bound, we 
will use a strong convergence result obtained in [12], combined with the idea 
that each uniformly integrable martingale Z7 essentially "counts" only the 
particles of corresponding velocity —7. 

2.3. The asymptotic shape and growth of the branching diffusion. The 
main result of this paper is the almost-sure rate of growth of particles which 
are in the vicinity of —jt in space and near K\ft in type position at large 
times t. For 7, re > 0, it can be shown that the limit 

(10) lim r 1 log E(N t (r, [«V*,oo))) 

t — >oo 



exists and takes the value 



A( 7 ,re):= , inf \E~ + A 7 - re 2 ^} 



Ae(A min ,o) 
(11) 

= p+ ^ - ^^/9{6 - 8r)(4adj 2 + a 2 {6 + K 2 ) 2 ). 

An outline of this expectation calculation is given in Section 3. Once again, 
we will find that the "almost sure" rate of growth of particles agrees with 
this "expected" rate exactly where there is growth in particle numbers. 

Theorem 3. Let 7, re > with A(7, re) 7^ 0. Under each P x > y law, the 
limit 

£>( 7 , re) := lim t" 1 \ogN t { T , [rev 7 *, 00)) 

exists almost surely and is given by 

A( 7 ,re), if A( 7 ,re)>0, 



(12) £>(7,«) = { 



-co, if A(7, re) < 0. 
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To prove the tricky lower bound of Theorem 3, which amounts to the 
major work of this paper, we will exhibit an explicit two-phase mechanism by 
which the branching diffusion can build up at least the required exponential 
number of particles near to — jt in space and K\ft in type position by large 
times t. 

During the first phase, over a large time t the process builds up an ini- 
tial excess of approximately exp(A(a)t) particles with spatial position at 
least —at, as is already known from Theorem 1. In this "ergodic" phase, 
"typical" particles found near —at in space will have drifted with a steady 
spatial speed of a whilst their type histories will have behaved roughly like 
OU processes with inward drift of \i\y for a certain optimal choice X(a) of 
parameter A. 

For the second phase, we will show that the probability any individual 
particle has at least one descendant that makes a "rapid ascent" in both 
space and type dimensions from initial position (0,0) to final position near 
(—(3t,K\/t) is approximately exp(— 6(/3, K)t), where 



and the time taken for this "rapid ascent" is an interval [0, r]. We show that 
this time r can be chosen such that 2[i\T ~ logt, and hence the additional 
time is asymptotically negligible in comparison with t. Intuitively, we will 
see that given an offspring that has successfully made such a difficult "rapid 
ascent," it will most likely have roughly had its type process behaving like an 
OU process with an outward drift of n\y and the Brownian motion driving 
its spatial motion will have had a drift A [corresponding to a real time spatial 
drift XA(y) that increases in strength as the type position y increases], for 
some optimal choice A(/3, k) of parameter A. The precise result required will 
be formulated rigorously as a large-deviation lower bound in Theorem 7 
of Section 5, and is proved using a "spine" change of measure technique 
intimately related to the Zt martingales. 

Combining these two phases and using independence of the particles, we 
can see that the number of particles near (— at,0) at time t that subse- 
quently proceed to have at least one descendant near is 
approximately Poisson with mean 



Optimizing for a fixed overall spatial speed 7, some calculus reveals that 



(13) 




y/0(0-8r)(a: 2 K*+4a0p 2 ) 
4a~6 



exp({A(a) -e(/3,/e)}t). 



(14) 




a+/3=7 
a,/3>0 
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with optimal parameters 

(i5) s=7 (^y and ~ p=i (jt^)- 

Thus we will be able to demonstrate an explicit two-phase mechanism pro- 
ducing the required number of particles, with this outline argument later 
guiding our rigorous proof. In addition, it is interesting to note that the 
optimal choices for A over each phase then also coincide at a single value 
\ = \(a) = \{p,K). 

An informative large deviation heuristic for the rapid ascent can also be 
found in Section 4, with this section also containing some essential optimal 
path calculations. We actually prove the two-phase mechanism for the lower 
bound of Theorem 3 in Section 5, although we defer proving the large- 
deviation lower bound until Section 7 after presenting the necessary "spine" 
background in Section 6. 

We prove the upper bound for the space-type growth rate in Section 9, 
again making crucial use of martingale results from Section 8. Similarly to 
the spatial growth case, we can find an upper bound using the Z+ mar- 
tingales, that is to say Nt(j; [«%/£, oo)) < exp{(£^ + A7 — K 2 ip^)t}Z^(t). 
However, as each martingale converges to zero, we must show that its 
exponential decay rate is (E^ — E^) before being able to optimize over the 
choice of A to obtain the required upper bound. 

Given Theorem 3, and noting symmetries, it becomes straightforward to 
retrieve the following: 



Corollary 4. For any Fcl 2 , define 

Nt(F):= 1 {{X u (t)/t,Y u (t)/Vt)&F}- 
ueN t 

If B C M 2 is any open set and CcK 2 is any closed set, then almost surely 
under any P x ' y 

liminf - logAf t (B) > sup D(j,n), 

limsup - logA/t(C) < sup D(j,k), 
t^oo t ( 7 , K )ec 

with the growth rate D(j,k) given at equation (12). 



We can also recover the almost sure asymptotic shape of the region occu- 
pied by the particles in the branching diffusion. 
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Corollary 5. Let BcR 2 be any open set. Almost surely, under each 
P x ^ law, 

f(n ro, if snB = 0, 

where 5cR 2 is the set given by 

5:={(7,k) GK 2 |A( 7 ,k) >0}. 

3. Some expectation calculations. This section discusses how the ex- 
pected growth rates given in the previous section may be obtained. For this, 
we use the "many-to-one" lemma (see, e.g., [8]) and one-particle changes of 
measure. In the process we shall start to gain valuable intuition into how 
particles within the branching diffusion behave, as well as seeing hints as to 
which are the "correct" martingales to use to prove the almost-sure growth 
rate results. 

For simplicity, we assume throughout this section that the branching dif- 
fusion starts with one particle at the origin in both space and type at time 
zero, unless otherwise stated. We also introduce a family of single particle 
probability measures P^a with associated expectations where, under 
Pu,A> V is an Ornstein-Uhlenbeck process with variance 8 and drift fj,, and 
£t = B(f A(rf s ) ds) where B is a Brownian motion with drift A. 

LEMMA 6 (Many-to-one). ///:1R 2 ^]R is Borel measurable then 

(16) E £ f(X u (t),Y u (t))=E e/2fi ^ex^j\(i h )dsy^ uVt )y 

Using the many-to-one lemma, and changing measure to alter the drift of 
Brownian motion, we see that 

E £ f(X u (t),Y u (t)) 

ueN t 

= E 0/2,o(exp^ R(ri s )ds)f(£t,r)t)) 
= E 0/2fi (V A «< exp (£{R{Vs) + y^(%)} ds 

= E 0/2)A ^ex P (-A6 + ^*{i?(r/ s ) + y^}^/fe ) 7 ?t )y 
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To perform a further change of measure on the OU process to get rid of the 
time integrals in the exponential of the expectation, we recall that 



and then 

E f(X u (t),Y u (t)) 

ueN t 

(17) =E fl/2iA (exp(-A& - ^r,t + # A • M^ m ) 

= E„ AiA (exp(-A& - VI Vt + E^t)m,Vt))- 
Note that the many-to-one lemma, combined with the branching prop- 
erty, immediately suggests how to get "additive" martingales for the branch- 
ing diffusion from single particle martingales — for example, taking f(x,y) = 
exp{\x + ip^y 2 } m equation (17) quickly leads to the martingale . 

We may now proceed to calculate the expected growth rates. However, 
for both clarity and brevity we will leave rigorous details to the interested 
reader, noting that the intuition we will gain from our rough calculations 
will later be invaluable in guiding our rigorous proof of the corresponding 
almost-sure growth rates. 

3.1. The expected rate of growth along spatial rays. We first give the 
outline of some calculations to find the rate of growth in the expected number 
of particles near — jt in space at time t. 

Using the formula from (17), for A G (A m i n ,0) and any e > we have 

E ZZ 1 {t- 1 X u (t)+ 1 d{-e,e)} 

ueN t 

< e^ +A ^- A£ )% AiA (e~^ 2 t ; y + 7 G {-e, e) 
> e (E-+x^+\e)t E ^ x L-^v 2 t ; y + 7 G (-£, e)\ 

where, with some abuse of notation that we shall continue to use throughout 
this section, we will abbreviate this to 

ueN t 

(18) 
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with the understanding that any subsequent arguments to identify expo- 
nential growth rates can readily be made rigorous by using the appropriate 
upper and lower bounds, and so on. 

Now, considering E~(X) := as a function of A, we have from (8) that 
A (7) = inf Ag ( Amhij o){-E _ (A) + A7} = £?~(A 7 ) + A 7 7, where A 7 satisfies 



/ n dE- ^ . , x (9-8r) 

(19) (A 7 ) = - 7 , hence A 7 = -71' 



OX y " " 1 1 V 9a 2 + Aa-f 2 

Of course, choosing this optimal A 7 value in (18) means that we must have si- 
multaneously maximized the expectation E (UAi ^(exp(— t/^ 7 /?)^ ~ — 7*)> and 
to confirm that this value is not exponentially decaying in t is now relatively 
straightforward. Under P^ Ai a 5 V is an Ornstein-Uhlenbeck process with an 
invariant measure given by the probability density, of the normal distri- 
bution N(0,9/(2nx)); and & = B( /* A(rj s ) ds), where B is a BM with drift 
A. Note also that by differentiating 

((Q e + (l/2)A 2 A + i?-i^K,^ = 

with respect to A, using self-adjointness, and observing that 4>\ oc (v^) 2 4>, 
we find that 



dE~ _(XAv^,v- 



X / A(y)(j) X (y)dy. 



Then almost surely under F > AtA ,A 5 

(20) 6 = But m *) & m ds ^x[A^ dy = x -^= d ^-, 

and so when we use the optimal A 7 value we get exactly the desired drift, 
since ^-(A 7 ) = —7. Then 



/ 

Jr 



'^ y2 ^(y)dy. 



In this way, we can obtain the exact rate of exponential growth for the 
expectation, 

lim r 1 log E(N t (r,^)) = Mi). 

The changes of measure used above are actually suggesting a great deal 
about the dominant particles that are found in the vicinity of a given ray 
in space. An alternative discussion of this expectation result, involving a 
dual approach via large deviation theory for occupation densities, can also 
be found in [13]. 
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3.2. The expected asymptotic shape. We give a rough outline of calcula- 
tions that will yield the correct exponential growth in the expected number 
of particles both near — jt in space and K\/i in type at large times t. Using 
the formula from (17) and abusing notation throughout in the same way as 
Section 3.1, we find that 

E Y 1 {X u (t) ~/t;Y u (t)>KVt} 

u£N t 

= ^,A(e"^-^^ + ^ i l {6 ._ 7t;??t > Kv/f} ) 
Now, from standard bounds on the tail of the normal distribution, 

P MA ,A(&~-7*;r/t>«v / t) 

(21) = P Ma ,a(^ > kv^)P Ma ,a(& ~ -lAm > «Vt) 

and, since tp^ + {pL\/0) = ^ , this yields 

E Y 1 {X u (t) ~ft;Yu(t)>K.Vt} 

ueN t 

(22) 



fix, 

• 2, 



Recalling that A(7, k) := mixe(X min fl){ E x + A7 — k ip^ }, simple calculus re- 
veals this infimum is attained at a A value of 



and using this optimal value in equation (22) will lead to the upper bound 

lim sup*" 1 logS 1 {x u (t)^- 1 t;Y u (t)> K Vt}^ A ^^y 
ueN t 

It is also clear from equation (22) that when minimizing E7 + A7 — , 
we simultaneously maximize the probability P^ Ai A(£t ~ ~lt\ r lt > Ky/t). In 
particular, to get a matching lower bound, we do not want this probability 
to have any exponential decay in time when we choose the optimal parameter 
for A. 

In fact, at least up to the exponential decay rate in time, it can be shown 
using large-deviations arguments that 



(6 ~ -7*; m > ~ ex P ( -^Vt ) . 
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Indeed, we immediately gain the required upper bound from (21). For the 
lower bound, consider the following heuristics where we break paths into two 
sections: normal ergodic behavior over large time period [0,t] followed by a 
rapid ascent out to type position K\Jt over a much shorter period [t, t + r] . 

(i) Ergodic behavior. Over a large time t, the occupation density of r/ 
will most likely have settled close to the invariant measure. Hence for large 
t, almost surely under Pa,/i a , 

(ii) Rapid ascent. Over a large time r, but where r = o(t), the probability 
that rj starts close to the origin and ends near to Ky/i, having followed close 
to the path y over the entire time period r, is roughly given by 



exp (~2^i) {y( s ^ + ^ xy ( s ^ 2 ds ) 



under the Fa,^ a law. See, for example [24], Chapter 6, or [4], Chapter 5.6. 
After some Euler-Lagrange optimization, the path 

r-sinh /j, x s 

y{s) = nVt—- 

sinn [j,\T 

gives Jq y(s) 2 ds ~ K 2 t/(2fix), with the probability of this path being roughly 
exp(-(fi x /9)K 2 t). 

Combining these two types of behavior, we can find paths with final positions 

rt+r / q K 2 \ 

T) t+T ~ K-Vt, Ct+r ~ Aa / rj 2 ds ~ Aa h - — t 

Jo \2^a 2^ A / 

and, moreover, when substituting the optimal A value of A(7, n) and sim- 
plifying, this actually gives £j+ T ~ —jt. Further, one of these paths occurs 
with a probability of roughly exjp(—(fix/6)n 2 t), and note that t + T ~ t since 
r = o(t). Thus we see that to exponential order, the probability P 



— jt;r]t > Ky/i) must be at least exp(— (iJ,^/9)K 2 t), as required. 

This heuristic argument can be made rigorous to prove, as claimed, that 

tonrtogtfl ]T l {Xll(t) <_ 7t; y u(i) > K V7 } =A( 7 ,«). 
\ueNt I 

If we scale all spatial coordinates by t~ l and all type coordinates by 
at time t, the expected asymptotic shape can be considered to be the 
region S := {(7, n) : A(7, k) > 0} where, on average, we have growth in the 
numbers of (scaled) particles. 
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4. Short climb large deviation heuristics. In this section, we give a heuris- 
tic calculation that suggests why the probability a single particle manages 
to have at least one descendant in the vicinity of (—/3t,K\/t) near time r 
is roughly exp(— G(/3, K)t) for very large t, where 0(/3, k) is given at equa- 
tion (13). For these heuristics, we will think of r as large and fixed, but of 
smaller order than t (later on, in our rigorous approach, we will choose r 
proportional to logi). We emphasize that the heuristics in this section are 
neither meant to be precise nor made rigorous, yet they will provide invalu- 
able intuition, guidance and motivation for our rigorous approach later on. 
Of particular importance will be the optimization problem that the heuris- 
tics suggest. Indeed, many of the exact calculations in Sections 4.2 and 4.3 
will be essential later in the paper. 

Suppose we start the branching diffusion with a single particle at (0,0). 
First, we wish to know the probability that there is at least one particle at 
time r that has a spatial position near —j3t having followed close to the path 
x{s) for < s < r and a type position near K\/i having closely followed the 
path y(s) for < s < r for t arbitrarily large. 

We recall from large deviation theory of Ventcel-Freidlin (see [24] , Chap- 
ter 6, or [4], Chapter 5.6) that the probability a single particle manages to 
follow closely both the type path y(s) and the spatial path x(s) for < s < r 
is roughly given by 

(24 ) e ^ /;( 9(s) + *,(»)) V^^ s ) 

when x(0) = 0,x(t) = —/3t, y(0) = 0, y(r) = K\/i and t is very large. This 
probability will typically be very small, but if such paths are followed by 
particles in the branching diffusion, we have to also take account of the 
large breeding rates that are found far from the type origin. 

If we let X(s) represent the numbers of particles in the branching diffusion 
that are alive at time s and have traveled "close" to the path (x(u),y(u)) 
for < u < s, then we can get a rough idea of how X might behave by 
considering the following birth-death process. 

4.1. A birth-death process. For given fixed paths x(-) and y(-), let M be 
a time-dependent birth-death process where at time s particles either give 
birth to single offspring with breeding rate A(s) given by 

X(s) = p + ry(s) 2 , 

or particles die with death rate /i(s) given by 

. . l/ n 9 . ,\ 2 1 x(s) 2 

^ = Ye{ y{s) + 2 y{s) ) + 2W 
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(Note that the probability the initial particle of this birth-death process sur- 
vives the entire time period [0, r] is consistent with the rough large deviation 
probability for the branching diffusion at equation (24).) 

An important quantity is the effective total death rate up to time t which 
is defined by u{s) := f${^{w) — \(w)}dw, so here 

v(s) = J(x,y,s) 

:= Jo{Te{ y{w) + 2 y{w) ) + 2^y~ ry{w) - p ) ds - 

The distribution for total number of offspring surviving, M(t), for the 
time-dependent birth-death process is well known, for example, see [14]. 
Then defining 

W T := e~ v ^ [\ + jT fi{s) e u{s) ds^j , 

U T :=l-e-^ T) W-\ 
V T :=l-W-\ 

we have 

F(M(t)=0) = U t , 

P(M(r) = n) = (1 - U T ){1 - V T )V?-\ n = 1, 2, . . . 

with EM(t) = e~ v ^ and E(M(t)\M(t) > 1) = W T . 
In our particular case, we have 

M(M(r)) = exp(-J(x,y,r)). 

Define the largest effective total death rate prior to time r by 

L(x,y,r):= sup J(x,y,s)>0. 

se[o,r] 

If we are in a case when L(x,y,r) is very large, suggesting a high chance of 
extinction, then 

(25) T{M{T) a 1] = i + T^w® is ~ Xt***-*™^- 

where K~ l := / Q T /i(s) exp(— {L(x, y, r) — J(x,y, s)}) ds. If there is at least 
one particle alive, we would then expect to have 

E(M(t)\M(t) > 1) ~ ^- 1 exp(L(x,y,r) - J(x,y,r)). 

Thus, we might guess that the probability any particles in the branching 
diffusion manage to make the difficult, rapid ascent along path (x, y) to finish 
up near (— /3t,Ky/i) can, very roughly, be estimated by exp(— L(x, y, r)). [To 
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help see this, try writing x(s) = tf(s) and y(s) = \/tg(s), thinking of /, g as 
fixed paths and recall that t is very large and r = o(t), then the role of K T 
in (25) is insignificant next to exp(— L(x, y,r)).] 

We might then further guess that the chance any particles manage to 
stay near position (—/3t, Ky/t) during a very small interval of time close to r 
should roughly look like 



exp(- inf L(x,y,r)j, 



where we permit all possible paths x and y satisfying x(0) = 0, x(t) = —j3t 
and ?/(0) = 0, y{r) = Ky/t for the fixed time r. (We will state and prove a 
precise lower bound that corresponds to this guess at Theorem 7.) 

4.2. Finding the optimal path and probability. We proceed to calculate 

'mIL(x, y,r) 
x,y 

over paths x and y satisfying x(0) = 0, x(t) = — fit and y(0) = 0,y(r) = K\/i 
for the fixed time r. 
We first note that 

(26) inf L(x, y,r) =inf sup J(x, y, s) > inf J(x, y, r) 
x,y x, y gg[0jT] x,y 

and we now proceed to calculate inf^y J(x, y, r). 

We can easily optimize over the choice of function x given y, finding that 

x(s) oc ay(s) 2 => x(s) = Xa / y(u) 2 du 

Jo 

where A is the constant of proportionality and must satisfy 

(27) X= a%yWd8> 
yielding 

1 r ±(s) 2 , I3 2 t 2 
■ as — 



1 r 

2 Jo 



ay(s) 2 2a^y{s) 2 ds' 



This is exactly as anticipated since, when following the path y in type space, 
the spatial position of a particle is following a Brownian motion with total 
amount of variance over period r given by a J Q T y(s) 2 ds. Hence, the proba- 
bility that a particle following the path y in type space will also be found 
near to fit in space at time r is roughly 

-I3H 2 

exp' 



2a Io y{s) 2 dsj' 
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Introducing the notation 



nd 

inf \l{y) + 



we are left to find 

(3 2 t 2 



2a So 2/O) 2 ds . 

(28) =infsup(l(y) - -aA 2 f y(sf ds - X(3t\ 

y x I 2 Jo J 

> sup inf (/(y) --aA 2 f y(sf ds - X(3t\ 
x y I 2 Jo J 

where the first equality is trivially true by maximizing the quadratic in 
A, the introduction of which conveniently removes the awkward integral in 
the denominator. Some further Euler-Lagrange optimization now gives the 
optimal path as 

, , . . ^sinh/i\s 

sinn /iAT 

where 



and then 



yJ9(9 - 8r - 4aA 2 ) 
Ma = : 



sup inf (/(y) - -aA 2 f y(sf ds - \pt\ 
x y I 2 Jo J 

sup|«; 2 tQ + ^coth/j, A r) - X(3t 



x 



The optimal parameter choice A (which depends on r as well as the model 
parameters) then satisfies 

(30) ^1 = kH ( £2^AI r ) = r ^ 

V ' aX V 2^ 2smh 2 H rJ Jo aV ' 

Then we have shown that 



2a Io Vx( s ) 2ds y I 2a Ioy( s ) 2ds 

= infsup{/(y) a\ 2 [ y(s) 2 ds — X/3t 

y x I 2 Jo 

> sup inf (/(y) - -aA 2 / T y(s) 2 - A/3tl 
A f I 2 Jo J 

>/( y ~)-I a A 2 f T y x (s) 2 ds - Xpt, 
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and, in fact, we see that the left- and right-hand sides of the above are equal 
by recalling (30). It follows that the preceding supremum and infimum can 
be freely interchanged, actually preserving equality at the inequality (28). 
Then, with the optimal spatial path 



(31) x\(s):=Xa / y\(u) du = —(3t 



flj sinh2/iAfi — 2fi\s 



sinh2^A'T" — 2/Uat ' 
and defining x := xs, y := y? , we have 
inf J(x } y,r) = J(x,y,r) 

= isup |k 2 Q + ^ coth^ A T^ - A/3 j - pr 

= 1 i" 2 (i + w coth ~ xp ) ~ pT - 

Finally, it is easy to check that J(x,y,T) = L(x,y,r), whence 

inf J(x,y,r) > infL(x,y,r), 

x,y x,y 

and, combining with equation (26), we have found that 
mlL(x,y,T) =MJ(x,y,r) = J(x,y,r). 

x,y x,y 

4.3. An important note on the optimal paths. As r — > oo, we have 
sup{/e 2 Q + y coth/i A r^ - A/3 J | sup{n 2 ^ - A/3} = k 2 ^ - A/3, 
where the optimizing parameters of the supremums also converge with 



Note the agreement with previous optimal values at equations (23) and (15). 
Then letting 

e(/3,re) := sup{rc 2 t/^ - A/3} 

(33) A 

_ k 2 v/^(0-8r)(a 2 K 4 + 4a^) 

and writing x := x^ and y ■= y\, we note that for all e,<5 > there exist 
f , p > such that for all t > and r > f 

exp I - inf J (x, y, r) J > exp(- J(x, y, r)) 

= exp f-t (k 2 (7 + 7^ coth Mx r J - A/3^ + pr 
> exp(-i(9(/3, /c) + e)). 
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Further (when k > 0), for all s £ [r — /U,r], 

y(s) > (k - 5)Vt, x(s) < -{p - 5)t. 

In particular, the paths stay close to the required positions for some fixed 
length of time with corresponding probability at least as large as required. 

5. Proof of Theorem 3. Lower bound. In this section we will state a pre- 
cise short climb probability result and show how to combine it with almost 
sure spatial (only) growth rates to prove the lower bound of the growth rate 
in Theorem 3. This will make rigorous the two-phase mechanism described 
in Section 2 and suggested by the expectation calculations in Section 3. 

The first phase requires knowledge of the almost-sure rates of growth of 
particles in the spatial dimension only. To this end, we will already make full 
use of Theorem 1 throughout this section, deferring its proof until Section 10. 

The second phase requires a lower bound for the probability that a single 
particle makes a rapid ascent in type-space over the time interval [0, r]. This 
is the lower bound found in the heuristics of Section 4, but we require some 
further notation before the precise result can be stated. Note, throughout 
this section, we will only be interested in the optimal parameter value A = A 
as introduced in Section 4.3. 

We wish to fix the relationship between sufficiently large t and r as 



(34) yJe/(2»- x )e^ = KVi 

and so define r = r(t) by 

r35) T(t) . = J(2^)- 1 log(2^/0), for 2fi- x t>6, 

y J w ' 1 0, otherwise. 

Recall the optimal paths (x,y) over s G [0, r], where 

. . . , r sinh lit s 

36 y(s)=KVt . — , 

sinh h x t 

(37) X (s)=aX [ S y( W )*dw = -pt Shlh2 ^ S - 2 ^ S 



sinh2/i^r — 2/x^r' 

with fixed end points y(r) = K\/t and x(t) = —fit. 
For large times t and 5, e > 0, let 

(38) A £ /(u):=\ sup \Y u {s) -y{s)\ < eVi; sup \X u (s) - x(s)\ < St \. 



se[0,r(t)] s6[0,r(t)] 

We will use the notation 

(39) Af:= U Af(u) 

u£N T(t) 
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for the event that there exists a particle in the branching diffusion that 
makes the short climb. Finally, recalling 0(/3, k) given at (33), we can now 
state the short climb theorem: 



Theorem 7. Fix any y\ > uq > 0, xel, and let So > 0. Then for any 
e, 5 > 0, there exists T > such that for all y £ [yoiyi]> 

for all t>T. 

We will prove Theorem 7 using a spine change of measure. This requires 
us to introduce the notation for the spine set-up in detail before proceeding, 
so this and further technical issues are postponed to Sections 6 and 7. 

Remark 8. We note that Theorem 7 is actually a stronger result than 
needed to prove Theorem 3 because we identify the specific paths followed 
by particles that are near position (fit, K\/t) at time t + r, rather than just 
considering the particle's positions close to time t + r. 

In combining the two phases, we will have a huge number of independent 
trials each with a small probability of success, intuitively giving rise to a 
Poisson approximation for a large number of successful particles. In fact, 
in our proof of the lower bound of Theorem 3 below, we will actually use 
the following result about the behavior of sequences of sums of independent 
Bernoulli random variables. 



Lemma 9. For each n, define the random variable B n := Ylu€F n ^-E n (u) 
where the events {E n (u) : u G F n } are independent. Let p n (u) := P(E n (u)) 
and S n := Y^u<=F n Pn(u) and suppose that, for some v G (1/2, 1), 



1 



(40) J2 ^ 

neN 



Then the sequence of [possibly dependent) random variables {B\, E>2, ■ . ■} 
has \B n — S n \ > {S n ) u for only finitely many n, almost surely. 

In particular, for any e > 0, there exists some (random) N £ N such that, 
with probability one, 

(41) > 1 - e for all n > N. 

Proof. For v e (1/2, 1), Chebyshev's inequality yields 



n) 
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and hence the Borel-Cantelli lemmas, combined with hypothesis (40), imply 
that 

\B n - s n \ > (s n y 

for only finitely many n, almost surely. Equation (41) now follows on division 
by S n , and noticing the assumption (40) implies that lmin^oo S n = oo. □ 

Proof of Theorem 3. Lower bound. Define / _1 (i) :=t — r(i), not- 
ing that both f{t)/t — > 1 and f~ l (t)/t — > 1 as t — ► oo. Also, for n G N and 
fx > 0, define T n := (n + We want to estimate the number of particles 
that are near the large position (—(a + f3)T n ,Ky/T^) during time interval 
[T n _i,T n ]. For this, we will consider particles that travel with a velocity 
—a over time period [0, f~ l (T n )\ before commencing their rapid ascent of 
(relatively short) duration r(T n ) to be in final position at time T n . Then 

inf N s ((a + (3-S)T n ;[( K - S)y/T^ )0 o)) 

s£[T n -i,T„] 

(42) > 22 1 {n se[T Tn] {Xu(s)<-(a+P-S)T n ;Y u (s)>(K-S)V^}} 

u£N Tn 

^ E 1 {Ni- K (u)>0} 
u£F£ 

where 

F% := {u G N f -, {Tn) :X u (r\T n )) < -oT n , Y tt (/- 1 (T n )) G [y ,yi]} 
and, for u G .F", 

iV^' K (u) := 1 {f| se[T _ Tn] {^(*)-^(/- 1 (T'„))<-(/3-5)T„;K ! ( s )>(« ; -<5)v / 7^}}- 

v£N Tn 
v>u 

We will now show that the sum at (42) grows as fast as anticipated: 

Lemma 10. For any e > 0, we may choose /i > such that there exists 
a random N G N where 

^ lo § E V^w>o}^ A ( Q )- (/ 3 ' K )- £ 

/or all n> N with probability one. 

Proof. We will be able to apply Lemma 9 given sufficient information 
about the growth of \F°\ and decay of the probabilities 

p^(u):=P(N^(u)>0\F r i {Tn) ), 
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where «£f*C Nf-i^y 

It follows easily from Theorem 1, f~ l {T n ) /T n — > 1 and the continuity of 
A (a) that 

loe|F a | £ 
S — 1 > A(a) - - 

for all sufficiently large n. 

The definition of N^' K (u) and spatial translation invariance implies that, 
for each u G F^, the rapid ascent probability Pn' K (u) depends only on the 
initial type position Y u {f~ l {T n )). 

For <5, fi > 0, define 

B S t '"(u):= p| {X u ( S )-X it (0)<-(/?-«5)t;y u ( S )>( K -5)Vt} 

se[r(t)-M,T(t)] 

and 

(43) B^:= |J B^(u). 

«e7V T(t) 

Recalling the comments of Section 4.3, there exist e',5' > and we may 
choose /x > sufficiently small, such that 

for all u G i 7 " whenever n is sufficiently large. Together with Theorem 7 and 
since l^(/ _1 (T n )) G [yo,yi] for it G F", this reveals 

logp^' K (?x) logpn(u) . nm s e 
^ ^ ^ ^ «J - - 

for all for u G i* 1 " and all sufficiently large n, almost surely. Then we may 
combine the observations above to obtain 

-Liog ]T pg'»>A(o)-8C9,«)-|. 

Taking this last line together the assertion of Lemma 9 at equation (41) 
gives the result. □ 

It is now straightforward to combine Lemma 10 with the inequality at 
(42) to see that, given e, 5 > 0, there exists /j, > and a random time T such 
that 

rhog N t ((a + (3- S)t; [(« - 5)Vt, oo)) > A(a) - 0(/3, k) - e 

for all t > T, almost surely. Since e and 5 can be taken arbitrarily small, 
using the optimal a and (3 according to equations (14)-(15), we find 

lim inf t ~ 1 log N t (7 , [kVI , 00 ) ) > A (7 , k) almost surely, 

t— >oo 
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as required. [It is also interesting to note that A = A s = A (7, k) from equa- 
tions (19), (23) and (32), so the optimal parameters are in agreement with 
those of the expectation calculations in Section 3 and the path large devia- 
tions in Section 4.] □ 

6. The "spine" setup and results. In this section, we describe how to 
construct an enriched branching diffusion with an identified "spine" or "back- 
bone" particle and discuss how to perform some extremely useful changes 
of measure (closely related to the additive martingales) that will essentially 
"force" the spine perform the short climb, whilst giving birth at an accel- 
erated rate to offspring that behave as if under the original measure. These 
spine techniques are at the very heart of our proof of Theorem 7 in Sec- 
tion 5. Spine ideas were first seen for branching Brownian motion in [3] and 
developed for Galton- Watson processes in [16, 18, 19]. Kyprianou [17] and 
Englander and Kyprianou [6], developed the technique for some families of 
branching diffusions; and more recently the spine approach has been signif- 
icantly improved in [8]. This approach uses several different nitrations on 
an enlarged probability space carrying the branching diffusion, and permits 
some very useful techniques and results to be developed. For example, "addi- 
tive" (many-particle) martingales can be represented as suitable conditional 
expectations of "spine" (single-particle) martingales and consequently there 
are clear interpretations for any changes of measure and all measures in- 
volved in our "spine" setup are probability measures with intuitive construc- 
tions. Following Hardy and Harris [8], we will first outline the notation and 
then describe the changes of measure. The notation described in this section 
is generalized to allow each particle u to have 1 + A u offspring, where each 
A u is an independent copy of a random variable with values in {0, 1,2,.. .}. 
The spine techniques developed in this paper could readily be generalized 
to such models. 

All probability measures are to be defined on the space T of marked 
Galton-Watson trees with spines; before defining precisely what this space 
is we need to set up some other notation. We recall the set of Ulam-Harris 
labels, Q, defined by Q := {0}UU Il6fJ (N) n 1 where N := {1, 2, 3, . . .}. For two 
words m,v6S1, uv denotes the concatenated word, where we take u0 = 0u = 
u. So 0, contains elements such as "0412," which represents "the individual 
being the 2nd child of the 1st child of the 4th child of the initial ancestor 
0." For labels u,v G SI the notation v < u means that v is an ancestor of u, 
and |u| denotes the length of u. 

We define a Galton-Watson tree to be a set rc!l such that: 

(i) G r, so there is the unique initial ancestor; 

(ii) if u, v G f2, then vu G r v G r, so r contains all of the ancestors of 
its nodes; 
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(iii) for all «Gt, there exists A u £ {0, 1,2,.. .} such that for j £ N, uj £ r 
if and only if 1 < j < 1 + A u . 

The set of all such trees is T, and we will use the symbol r for a particular 
tree. As our work concerns branching diffusions we shall often refer to the 
labels of r as particles. Note that for the binary branching mechanism in this 
paper, P(A U = 1) = 1; of course, here there is only one r £ T — the binary 
tree. 

A Galton- Watson tree by itself only records the family structure of the 
individuals, so to each individual u £ r we give a mark (X u ,Y u ,o~ u ) which 
contains the following information: 

• cr u £ [0, oo) is the lifetime of particle u, which also determines the fission 
time of the particle as S u := J2v<u a v- We may also refer to the S u as 
death times; 

• the function X u (t) : [S u — a u , S u ) — > M. describes the particle's spatial mo- 
tion in K during its lifetime; 

• the function Y u (t) : [S u — cr u ,S u ) — > E describes the evolution of the parti- 
cle's type in R during its lifetime. 

For clarity we must decide whether or not a particle is in existence at its 
death time: our convention will be that a particle dies "infmitesimaHy be- 
fore" its death time — this is why X u and Y u are defined on [S u — a u , S u ) and 
not [S u — cr u ,S u ] — so that at time S u the particle u has disappeared and has 
been replaced by its two children. 

We denote a particular marked tree by (r,X,Y,a), or the abbreviation 
(r, M), and the set of all marked Galton- Watson trees by T. For each 
(r,X,Y,a) £ T, the set of particles alive at time t is defined as N t := {u £ 
t : S u — a u < t < S u } . For any given marked tree (r, M) £ T we can distin- 
guish individual lines of descent from the initial ancestor: 0,u±,U2,U3, ... £ 
r, where Uj is a child of for all i £ {2,3, . . .} and u\ is a child of the 
initial individual 0. We call such a line of descent a spine and denote it by 
£. In a slight abuse of notation we refer to £j as the unique node in £ that is 
alive at time t, and also for the position of the particle that makes up the 
spine at time t; that is, £j := X u (t), where u £ £n Nt- However, although the 
interpretation of should always be clear from the context, we introduce 
the following notation for use where some ambiguity may still arise: 

• nodet ((t, M, £)):=« if u££ is the node in the spine alive at time t. 

It is natural to think of the spine as a single diffusing particle £t, or, strictly 
speaking, the pair (£t,7?t), where r\ t is the type of the spine at time t. 

We define nt to be a counting function that tells us which generation of 
the spine is currently alive, or equivalently the number of fission times there 
have been on the spine: 

n t = |node t (£)|. 
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The collection of all marked trees with a distinguished spine is the space 
T on which our probability measures will eventually be defined, but first we 
define four nitrations on this space that contain different levels of information 
about the branching diffusion. 

• Filtration (T)t>o- We define a filtration of T made up of the cr-algebras 

T t := a{{u, X u ,Y u ,a u ) : S u <t; 

(u,X u (s),Y u (s):s G [S u - a u ,t]) :t € [S u - a u , S u )), 

which means that T is generated by the information concerning all par- 
ticles that have lived and died before time t, and also those that are 
still alive at time t. Each of these cr-algebras is a subset of the limit 
Too := a(\J t > T t ). 

• Filtration (T)t>o- We define the filtration (T)t>o by augmenting the 
filtration T with the knowledge of which node is the spine at time t; that 
is, (T)t>o := n odet(£)) and T^ := cr(\J t>0 Tt), so that this filtration 
knows everything about the branching diffusion and everything about the 
spine. 

• Filtration (Gt)t>o- (Gt)t>o is a filtration of T defined by Gt ■= cr(£ s : < 
s<t), and Qoo := a([j t>0 Qt). These a-algebras are generated only by the 
spine's motion and so do not contain the information about which nodes 
of the tree r make up the spine. 

• Filtration (Gt)t>o- As we did in going from Tt to T we create (Gt)t>o 
from (Gt)t>o by including knowledge of which nodes make up the spine: 
(Gt)t>o ■= f(&,node t (£)) and := cr(\J t > Q Gt)- This means that Gt also 
knows when the fission times on the spine occurred, whereas Gt does not. 

Now that we have defined the underlying space and nitrations, we can de- 
fine the probability measures of interest. We let the typed branching diffusion 
be as described in Section 1.1, with the probability measures {P x ' y : x, y E M} 
on (T,Too) representing the law of this typed branching diffusion when ini- 
tially started with a single particle at (x,y). 

We recall from [18] that, if / is an J^-measurable function, we can write 

(44) /=E 

ueN t 

where f u is ^-measurable. Now we can extend P x ' y to a measure P x ' v on 
(T,Too) by choosing the particle that continues the spine uniformly each 
time there is a birth on the spine; more precisely, for any / G mTt with 
representation like (44), we have: 
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We construct the .^-measurable martingale £(t) as 

£( t ) := v + {r]t)e f*{R(Vs)+l/2X*A( Vs )}ds-E+t x TH - f*R(Vs 
(45) ><e A€t-l/2A»/ (J *A(^) < fa 



V 



Observe that this is a product of single-particle martingales, details of which 
can be found in [17] or [10]. One can think of these as /i-transforms of the P- 
law of the spine: the first makes rj an outward-drifting Ornstein-Uhlenbeck 
process with drift parameter the second increases the breeding rate on 
the spine to 2R(-); and the third adds a spatial drift to £. 

Using the martingale ((t) we may define a measure Q^' y on (T,.^) by 



(46) 



dP x <y 



C(0) v+(y) 



And since ((t) is a product of /i-transforms, under Q^' v the process may be 
re-constructed path-wise according to the following description: 

• starting from spatial position x and type y the spine (£,t,i]t) diffuses spa- 
tially as a Brownian motion with infinitesimal variance A(rjt) and infinites- 
imal drift XA(r]t); 

• the type of the spine, rjt, begins at y and moves in type space as an 
outward- drifting Ornstein-Uhlenbeck process with generator 

6 d 2 d 

• the spine branches at rate 2R(r]t), producing 2 particles; 

• one of these particles is selected uniformly at random; 

• the chosen offspring repeats stochastically the behavior of its parent; 

• the other offspring particle initiates a P V -BBM from its birth position 
and type. 

The change of measure (46) projects onto the sub-algebra Tt as a condi- 
tional expectation: 



dP x >y 



-Xx 

P^( v +(r, t )2^- E ^ t ), 



and it is a short calculation using the methods of, for example, Hardy and 
Harris [10] to show that: 
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Theorem 11. If we define ^ A 
J-oo that satisfies 



^ V \f<x>> then Q\ is a measure on 



17 



dP x >y 



zUt) , 



Moreover under Q x x ' y , the path-wise construction of the branching diffusion 



is the same as under 



Although the path-wise construction of the branching diffusion is the same 
under Q^' y and Q^' y , only the measure Q^' v "knows" about the spine. It is 
clear, however, that we have Qf^ y (A) = Q^' V (A) for any A G J 7 ^. 

Under the measure Q^' v only the behavior of the spine is altered, and 
combining this observation with conditioning on the spine's path and fission- 
times gives us a very useful representation for (t) under Q^' v that we shall 
refer to as the spine decomposition: 

(47) Q^Z+WI&o) = £ v+insJtX**-** 8 " + vt(vt)e^- E tK 

u<£t 

Throughout the rest of this article we will refer to the two pieces of this de- 
composition as the "sum term" and the "spine term." This decomposition 
is discussed in detail for a wide variety of branching diffusions in [9], but 
to derive it we simply note that the contributions to Zf{t) from the sub- 
trees that branch off the spine have constant Q^-expectation because they 
behave as if under the original measure P, and we know that Zt(t) is a P- 
martingale. The spine decomposition reduces many calculations about the 
behavior of Z^(t) under Q^ y to one-particle calculations about the spine, 
and this observation is exploited in the spine proofs of £ p -bounds for some 
families of additive martingales in [9]. 



7. Proof of Theorem 7. The short climb probability. With the spine 
foundations firmly established in Section 6, we may proceed with the proof 
of the short climb probability lower bound from Theorem 7. 

First, recall definitions (38) and (39), where A e t ' S is the event that there 
exists a particle that makes the short climb along optimal path (x,y), and 
Al' S (£) is the event that the spine makes the short climb. Note that e controls 
the proximity to x and 5 the proximity to y. Importantly, we will only be 
interested in taking A = A throughout this section, although we will usually 
just write A for notational simplicity. Also recall throughout that t and r 
are related through (0/(2n\)) exp(2/i^r) = K 2 t. 



A TYPED BRANCHING DIFFUSION 



29 



Proof OF Theorem 7. The key step in the proof of this is the following 
use of the spine change of measure: for any function g : R + — > ]R + we have 

(48) > Q^ffef; sup Zt(s)<g(r)) 
ygirr^UfiO; ^p Zt(s)<g(r)). 

V se[0,r] / 

Essentially we just have to make the "correct" choice for both A and g in 
expression (48), although there will still remain a number of technicalities 
to resolve. 

The first idea is to ensure the (originally rare) event A t ' actually occurs 
under the new measure (Qf^' y by making the spine follow close to the required 
path (x,y); this is achieved by choosing the optimal value A for A and choos- 
ing t to be on the natural time scale it would take the spine to reach position 
K\ft. In particular, this choice will mean that in the first line of the above set 
of inequalities there is no significant loss of mass when replacing the event 
A £ t ' 5 with Af S (£). Next, we wish to choose the smallest possible g that will 
still leave some positive probability on the last line of the above argument. 
Hence, we wish to identify the rate of growth of the martingale under 
Q?' 27 , and this will essentially be governed by the contribution from the spine 
itself. 

With this is mind, and recalling the various properties of the optimal 
paths and parameters from Section 4, for eo > we define 

*.M := «p((tf + ^ + f ) (^> 2 - - W> A V + A,)) 

and recall from (35) that the scaling between t and r is fixed throughout, 
where K 2 t = (#/(2/^))e 2 ^ r for large t, hence t + r ~ t. Note that since we 
are only considering the optimal value A = A, we have 

Then from (48) we have 

(49) P^iA^^g^rr'Q^UfiOl sup Z+ (a) < g £o (r)Y 

V sS[0,r] / 

Our strategy for the rest of this proof is to show that the Q^ ,s/ -probability in 
(49) is at least some e' > for all sufficiently large t, uniformly for y E [yo, yi], 
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so that the decay rate part of (49) matches the desired rate in the statement 
of the theorem. 

Conditioning on the spine's path and birth times, Goo, and then making 
use of some standard properties of conditional expectation we have 



jx,y( A e,8 



A t ' (0; SU P z x ( s ) <9e {r] 

s6[0,t] 



frMQ^K-'CO; ^p Z+(s)<g E0 (r) 



vf,y 



se[o,- 

sup Z+(s)<g £0 (T] 

s£[0,r] 



since Af S (£) is (?oo-measurable. We next observe that, conditional on Goo, 
we can write Z^(t) as 

(50) Zt(t) = e -Wt*+*> ( £ e- E P»Z^\t - S u ) + /(*)) , 



r(«) 



where the are independent copies of started from a single particle 
at (£s u ,VSu)'i and /(t) is the contribution to Z~^{t) from the spine, which, 
conditional on Qoo-j is & known function of t. Now if we could show, for 
< e < e , 

sup f(s) < ^ and sup (Z+(s) - /(*)) < ^ll, 

se[0,r] z se[0,r] * 

where /(t) := y 2+Xx ^f(t) , we would have sup sG [ r ] ^(s) < <fe (r). Hence, 
defining Z^{s) : = Z^(s) — f(s), we have 



.r.lj 



sup Z+(s)<g £0 {r) 
>e[o,r] 



> 



?/ \ ^ 9eo\ T ) \ ^ 9eo\ T ) 

sup f{s)<— ^ — ; sup -Z^(s)< 



se[0,r] 



se[0,r] 



,1 



A I X' (0 W, 6 [0,r] /(«)<J«oW/2} V * 



sup -Z^(s) < 

re[0,r] 



2 



since, conditional on C/oo, the supremum of / on [0,r] is known. 

We see from (50) that, conditional on Goo, Zt(t) is a submartingale. This 
is because the Q^' y -conditional expectation of each of the z| in the sum 



(51) 



e -(^V+A*) J2 e- E t s -Z^\t - S u ) 
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9e ( T ) 



is constant, so the expectation of the sum cannot decrease, and in fact this 
expectation increases every time there is a birth on the spine. Then by 
Doob's submartingale inequality we have 

sup Zt(s)< 9 -^ 

se[0,r] L 

= 1-Q^f sup Z+{s)> 

Vse[0,r] 

We must note here that the expectation on the above line is not a priori 
finite. However, the expectation of each term in the sum (51) is bounded by 
sup sg [ 0jT ] f(s), which we have control over via an indicator function and so 
we do not have to worry about this expectation blowing up. 

So we need to show that for all sufficiently large r and all y £ [yo, yi], 

or(^M ran{raP . tM / w < 9(oW/2 ,(i- i ^«r(4 + wiff»)))>^ 

and hence also 



Afit); sup Z+(s)<g E0 (r)) >e' 



s6[0,' 

as required. This will follow by combining both parts of the following result. 



Lemma 12. Fix y\ > y > and e > i > 0. 

(i) For all sufficiently small e, 5 > 0, there exists some e' > and T > 
such that for all y S [yo, yi] and all t>T, 

Q?»(V(0; sup f( s )< 9 J^l) >£ '. 

V s6[0,t] 1 / 

(ii) As t — > oo, 

®l' y (^®7(ZUr)\Gooy, sup f\ s )< 9 lM\^0 
\9e {T) s e[o,r] 2 / 

uniformly over y S [yo,yi]- 

Then we have shown that, for any eo > 0, yi > yo > 0, and sufficiently 
small £, (5 > 0, there exists a T > such that, for all y G [yo, yi] and a h t>T, 

t- 1 iogP x 'y(Af)>-(e(p,K)+e ). 
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Finally, we observe that the probability P x ' y {A £ t ' 5 ) is trivially monotone 
increasing in both e and 6, and so it follows that if the result is true for all 
sufficiently small e and 5, it is in fact true for all e, 6 > 0. This completes 
the proof of Theorem 7. □ 



Proof of Lemma 12 (i) . We will prove Lemma 12(i) in a sequence of 
other lemmas, using a convenient coupling for the spine's type process. 
First recall that, under Q^' v , r] s solves the SDE 

drj s = V9dB s + nxiy ds, 

where B s is a Q^-Brownian motion. Noting that d(e~^ xS r] s ) = e~^ xS V9 dB s , 
we can construct e~^ xS r] s as a time-change of a Brownian motion with 



e 



-^ s Vs -r, = V9 f dB w = J JLb(1 - e~ 2 ^ s ) 

JO V 2 /"A 



where B is also a Q^' y -Brownian motion started at the origin. 

In this way, we will construct rf under P from a Brownian motion B y 
started at y^j2i7\[9 where, for s G [0, oo), 



rfU) = J— e ^ s By(l - e~ 2 ^ s ). 
V 2 A*A 

To construct simultaneously all type processes rf under the same measure 
P, we first construct the process B yo as an independent Brownian motion 
started at yo\/2fi\/9. Second, we construct the process B Vl by running an 
independent Brownian motion started at y\\/2ii\/6 until it first hits the 
path of B Vo , at which point we couple the two processes together. Next, for 
any other y G (yo,yi), we run an independent Brownian motion B y until it 
first meets with either the process B yo below or B yi above, at which point 
we couple it to the process it first hits. 

Finally, we construct all the corresponding spatial processes under P 
from a single Brownian motion W by defining 



s 



(52) £ y {s) = W[aj r, y {wfdwj+\a J rj y {w) 2 dw, 

where W is started at x and is independent of the B y processes. 

Constructed in this way, for each y G [yojyi]j the P-law of (£, y ,r] y ) is the 
same as the Q^' y -law of (£,r)). 

Fixing fj, G (0, 1) and K > maxjyi, 1}, we define the events and stopping 
times 



A y := \ B y (s) G 



e -, £ 

1 ,H 

2k 2k 



>e(l-M] 



T := inf{t : B yo {t) = 0}, T K := inf {t : B yi (t) = K}, 
A £>K := Af n Af n {T > 1} n {T K > 1}. 
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Then, clearly P(A £t x) > and, on the event A £t K, the coupling gives 



0<rf°(s) <rf(s) <T] yi (s) <KJ e^ s , 

for all s > and y £ [yo,yi\- Note that our construction also ensures that if 
event A yo n A yi occurs then so must A\ for any y £ [yo, yi] , hence A\ D A £j k- 
□ 



Lemma 13. Let e > 0. On event t/iere exists a deterministic time 

sq = s (e) > such that for all r > s , 



sup \if(s) — y(s)\ < eVt, 

sS[0,t] 



for all y £ [yo,yi]. 



Proof. Set si = — T^Tog/x and then, on event A e ,k, f° r all r > s > si 
we have 



- 2k V 2 fix 



for all y £ [yo,yi\- Writing 



9 



5 W = fc e "\i- e — 



< 



we see that there exists S2 = S2(e) > such that, for r > s > S2, 

e / 



2^A 



< 



P^A« 



2k y 2^ 



Taking S3(e) = max{si, S2(e)} now yields 

e / 



(53) 



\r] y (s)-y(s)\ < 



pMa. 



K V 2 MA 



<ey/i 



for all r > s > S3 and all y £ [yo, y±]. 

Now consider s € [0, S3]. On A £i k we have 



W{s)-y{s)\<^—e^(l + K), 

and hence for some 84(e) > we have \rj v (s) — y(s)\ < E\ft for all r > S4, all 
s£ [0,53], and all y£ [yo,yi]- Taking sq(s) =max{s3,S4} yields the result. 
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Lemma 14. Let 5 > 0. Then for all sufficiently small e, there exists a 
deterministic tq = To(e,6) > such that, on A £ ^k, we have 



(54) 



sup 

sS[0,t] 



rf{w) dw — / y(w) dw 



< 5t 



for all t > tq and all y £ [yo,yi\- 



Proof. Given any 5 > 0, we first fix an e > sufficiently small such 
that e(2 + f ) < f I this yields a corresponding S3 = 53(e), which is chosen 
as at equation (53). Given this S3, we find n = ri(e, 5) > such that, for all 

r > n, 



(K 2 + l) r — e 2 ^ w dw < -t. 
Jo 2/i A 4 



We now set To = To(e,5) = max{s3,Ti}. With this choice of e and To, we 
proceed to show that the inequality (54) is satisfied. Note that tq is deter- 
ministic and independent of y. 

From equation (53) we see that, on A £ ,k and for s > S3, 



rS3 rs / £ 9 

i] y (w) 2 dw> / rf{w) 2 dw + / I y(w) J 1 

Jo J s 3 \ K y 2/j.x 



dw 



> / y(w) dw 



y(w) 2 dw — 2 



e9 

1 

3 2k Ma 



dw 



rs rss a K f 

> / y(w) 2 dw- / — e 2 ^ w dw-{2e) — 
Jo " Jo 2\x\ 2nx 

rs J 

> / y(w) 2 dw--t 
Jo 2 



for all r > tq and all y G [yo>2/i]- Similarly 



(\ v (wfdw< ( Si rf{w) 2 dw+ f (y(w) + -\hr~e^ w 

Jo Jo Js 3 \ « V 



(fu> 



dw 



+ 



2 + 



2^ A 



< / y{w) 2 dw + K 2 



< / y(w) 2 dw+ -t 



3 ^e 2 ^dw + e(2 + £ -)^ 
2fi x \ kJ 2/i A 
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for all r > tq and all y G [yo,yi\- Finally, for s G [0, S3], on A Et K we have 



rf(w) 2 dw — j y(w) 2 dw 



< / r) y (w) 2 dw+ / y(w) 2 dw 



< (K 2 + l) / — e 2 ^ w dw<5t 
Jo 2^i A 

for all r > ro and all y G [yo 5 Hi}- D 

Lemma 15. Ze£ <5 > 0. T/ien /or all sufficiently small e > 0, there exists 
F-almost everywhere on A Ei k cl random time Sq = Sq(5,e) < 00 such that 



sup 

ss[o,t 



£ y (s)-Xa I rf{w) 2 dw 
for all y G [yoj^i] o,nd all r > Sq. 



<St, 



Proof. Given 5 > 0, choose any S',S" > such that 5'(\f3/\\ + 5") < 5. 
Recalling the construction of £ y at (52), we see from standard properties of 
Brownian motion that there almost surely exists some S\ = Si(<5') < 00 such 
that 



- sup \W(s)\<5' for all f >S X . 

* se[o,t] 



Then 

(55) 



sup 

«6[0,t] 



W\a r] y (w) 2 dw 



<6'la rjy(w) 2 dw 



for all r such that a JJ" n y {w) 2 dw > Si, and by the coupling construction, 
on A £j k this is true for all y G [yo,yi] if a J^rj yo (w) 2 dw > Si. Then there 
exists (P-almost everywhere on A £i k) a random time S2 = 5*2 (5 ; ) < 00, which 
depends on B Vo and Si, such that a J Q r rf{w) 2 dw > Si for all y G [yo, yi] when 
r>S 2 . 

Now by Lemma 14, given 5" and a sufficiently small e, there exists a 
deterministic ro = To(e, 5") > such that, on A e< K, 



(56) 



a / r\ y (wydw<a\ y(s) 2 ds + 5"t 



X 



+ 5"\t 



for all r > ro and all y G [yo,yi\- Combining the inequalities at (55) and (56), 
we now see that, for t > Sq = So(e, 6' , 5") = max{S2,ro}, 



sup 


^(s)-Aa f'rf(w) 2 dw 


= sup 


W^a J* n y (w) 2 dw^j 


< St 


S[0,t] 


Jo 


se[o,r] 







for all y G [yo,yi\. □ 
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On combining Lemmas 14 and 15 and recalling the definition of optimal 
path x at (37), we obtain the following: 

Lemma 16. Let 5 > 0. Then for all sufficiently small e > 0, there exists 
P-almost everywhere on A Et x a random time Sq = So(5,e) < oo such that 

sup \£ y (s)-x(s)\<St, 

se[0,r] 

for all y € [yo,Ui], and all t > Sq. 

We may now draw everything together to finish the proof of Lemma 12 (i). 
First we observe that since A < 0, on event (£), 

sup e -«?/ 2 +^) exp(^+^ + - E+s) 

s6[0,t] 

< e -W+y 2 +Ax) exp (,/ ; +( K + e f t + \(-l3 - 8)t), 
and so, given £q, we can choose first 5 and then e sufficiently small so that 

Af(£)c{ sup fc)<%P) 
Ue[o,r] 1 J 

and, from Lemmas 13 and 16, there exists a random time T = T(5,e) < oo 
such that on A Et x we have 

sup \rj y (s) — y(s)\ < eVi 

s6[0,t] 

and 

sup |£ y (s)-x(s)| <<ft 
se[0,r] 

for all r > T and all y G [yoj That is, A E) jf D {T < r} C Af ,<5 (£ y ) for each 
V ^ [yoiZ/i]) with the slight abuse of notation that 

Af(e) = \ sup \rjy( s )-y( s )\<eVi; sup |^( S ) - x(s)\ < St). 

Lse[0,r(t)] sG[0,r(t)] J 

Note also that P(t4 £i x) > e' for some e' > 0. 

Combining the above, for any y £ [yo,yi] we have 

q^(V(6; sup /( S )<^M)=Q^(^ a (a)=p(^ a (a) 

V se[0,r] 2 / 

>P(I £ , A ';T<r)^P(i £ ^) 

as r — ► oo, as required. □ 
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Proof of Lemma 12 (ii) . Consider the expectation of the "sum term." 
We have 



S'Va 



Hence 



(i+(r)|g 0O ) = e-W^ 2 +^)Q^ £ e- E t^Z^(t-S u ) 

< e -(^y 2 +\z) nrmax{e iPtv(S u f+xaSu)-E+S u . u< ^} 

<n T sup /(s). 

se[0,r] 



K' y {^&x y (Zt(T)\Goo); sup /(.)< ^ 



3eol r J se[0,r] 
(57) < Q^fn T ^4; sup /(,) < ^ 

< e^-^Q^K), 

and we can now calculate Q^' v (n T ) = Q^'^Q^'^fvlC/oo)), where Goc the 
cr-algebra generated by the path of the spine (not including the birth times). 
Conditional on G^, n T is a Poisson random variable with mean given by 
/ r 2(ri] 2 + p) ds, and using Fubini's theorem we have 



f x ' y (^2(r V 2 + p)ds 



2rQ$ y (r}l)ds + 2pr 
— — + y 2 e 2 ^- — + y 2 ) + 2pr 

—K 2 t+ y / A t + o{r). 
/"A V 



So the Q^'^-expectation of n T only grows linearly in t. Then since £q — £q > 0, 
the expression at (57) tends to as t — > oo. Moreover, the expectation at (57) 
is bounded by the Q^' V1 -expectation, and hence the convergence is uniform 
over y € [yo,yi], as claimed. □ 



8. Martingale results. In this section we recall some existing and prove 
some new martingale results that are intermediate steps in the proofs of 
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Theorem 1 and the upper bound of Theorem 3. We recall from [13] that E\ 
[also written E~(\)] and A(7) are Legendre conjugates with 

(58) A( 7 ) = inf IS- (A) + A 7 }, BT (A) = sup{ A( 7 ) - 7 A}. 

A<0 7 >o 

If, for A m i n < A < 0, we write 7 a for the 7 value which achieves the supremum 
on the right-hand side of equation (58), then the functions A i— > 7 ^ from 
(— A m i n ,0) to (0,oo), and 7 i— > A 7 from (0,oo) to (— A m i n ,0) are inverses of 
each other and, of course, A 7 is the A value which achieves the infimum on 
the left-hand side of equation (58). In addition, we note that 



(59) 7, = 



d\ v ' V - 8r - 4 « A2 ' 
that E~(\) and A( 7 ) are convex functions, and that 

c(9) = sup{ 7 : A( 7 ) > 0} = inf {-ET (A) /A : A min < A < 0} 



(60) 
where 
(61) 



inf{c A : A min < A < 0} = c- , 



-E^ /A and 



m - "V a(e + A P y (Amin ' 0) - 



A formula for c(9) is given in equation (9). The following fundamental con- 
vergence result for the martingale was first partly proved in [13], but 
also see [9] for a more complete proof using "spine" techniques. 

Theorem 17. Suppose AG (A m i n ,0]. 

(i) // A 6 (A(#),0] ; the martingale is uniformly integrable and has 
an almost sure strictly positive limit. 

(ii) If A < \{9), then Z~^ (oo) = almost surely. 

The following convergence result was proved in [12] using martingales 
based on Her mite polynomials. 

Theorem 18. Let A e (A(0),O] and a < 1/4. For each P x ' y starting law 
and every continuous bounded function / :Mnl, we have 

/( y u( t ))e ^(^ + A(^(t) +c - t ) _>/ Z A -(oo), 
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where 

(62) /„:= (^) V4 Jj{y)e a y 2 e^y 2 4 ) {y)dy 

and <j>(y) is the standard normal density. 

In this paper, we require a corollary to this theorem which specifies more 
precisely which particles contribute to the final limit. 

Corollary 19. Let A G (A(0),O] and a < 1/4. For each P x > y starting 
law and every continuous bounded function /:Ri— ► R, we have for every 
e>0 

(63) £ /(F u (t)) e ^(*) 2 +^(*)- E A*i { | Xu(t)/mA | <E}i:r /oZ A -(oo) 

where 7a = — -^E~ (X) and /o is gu>en ai equation (62). 

This last result will enable us to show in Section 10 that the almost sure 
growth rate is at least as large as the expected growth rate, .0(7) > A(7). It 
is easy to see from Corollary 19 that when Z^(oo) > 0, there must exist at 
least one particle near to —j\t in space. Further, because of the decay rate 
of each term in the sum over particles at equation (63), it will be relatively 
straightforward to improve this to get the required exponential numbers of 
particles, exp(A(7)t), near — 7a t for large times [as long as Z^(oo) > 0]. 

The following result concerns the rate at which the martingales and 
Z7 converge to zero. 



Theorem 20. Let A e (A m j n ,0). For every starting law, P x > y ', 

^(cf-ct) a.s. 



log^(t) 
where is given at (5), and 



f 



5(9), if A min <A<A(0), 

c-, if ~x(e)<x<o. 



Corollary 21. If \€ (A min ,0), then Z^(t) — ► P x > y -almost surely. 

The rate of convergence of the Z^ martingale in part (i) of Theorem 20 
is crucial in Section 9 to obtain the upper bound on the almost sure growth 
rate, D(^,k) < A(7, k). We also comment that if Corollary 19 were true 
for all a < V>a > then we could have gained this upper bound at that point. 
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Although Corollary 19 is only proven for a < 1/4 (where we can utilize 



suitable Hermite expansions), we conjecture that it holds for all a < ifi~f. 



Proof of Corollary 19. Let e > be small, n ■— X— e, A,/i e (A(0),O), 
f be a positive, continuous bounded function, a < 1/4 and note that 7^ > 7a- 
Then we have 

E /(^(*))^ ( ^ + ^ (t) -^V.w<-7,t } 

< e(fi --^-- £7M )t ^ /(yu(t))e ay tt (t) 2+ ^( t) - E -* 1{Xu(t)< _^ } 



ay„(t) 2 +M^n(t)-S-^ e -(£;- -S,7+(A-m)7m)«. 



< E /(y«(t))e 

Recall that E~(X) is convex with -j^E~{X) > and ^^"(A) = 7a, so, from 
the Taylor expansion, 



dX 

(m-a) 2 d 2 



ox 2 



E-{X) + o{{n-XY). 



Then taking e > small enough so that E^ — E^ + (A — /u)7^ > 0, and using 
Theorem 18, we find that for any 5 > 

limsup £ /(y u (t))e^W 2+A ^W-^*l {Xu(i)< _ (7A+5) , } =0. 

u&Nt 

Similarly, we can show 

limsup £ /(y(t))e^W 2+ ^W-^ 4 l {x „ w> _ {7A _ 5)4} =0, 

and hence the only contribution to the limit comes from the particles near 
— 7a i in space. Combining this with Theorem 18 we have 



aY u (tf+\X u (t)-E~t 

t— >oo 

ueN t 



/ ^(oo)= t hm J2 f(Xu(t))e 

u&N t 

= lim £ /(F u (t))e^W 2 +^W-^ t l {|x4i)/t+7A|<5} . 

«G7Vt U 

Proof of Theorem 20. We use a useful technique brought to our 
attention in [22]. Let p£ (0,1) so that, by Jensen's inequality, Z^(t) p is a 
supermartingale; then for u,w>flwe have 

(u + vf <u p + v p , 
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and hence 



e V>^(*) 2 +Apf u (t)+c±t) 



ueN t 



< [ V e^* y «(*) 2 +pA(X„(t)+c- A t) j gPA(c 



M£N t 

For any e > 0, Doob's supermartingale inequality says 

EZ£(s)P 



P{ sup Z±(w) p >e p )< 



s<w<s+t 



s p 



V ueN s ) 



and then 



P( sup e dw Z±(w)>e 

s<w<s+t 



<P[ sup Z±(w) p >e" p<5(s+ V 

\s<w<s+t 



c x~ c P x)+ 5 ) s 



<: £ -p e pSt( E e P^x Y v(. s ) 2 +pHXu(s)+ CpX s) \ e p(H 
V ueN a / 

Now, if we can choose p £ (0, 1) such that A(c^ — c~ A ) + 5 < and pV'a < - 
^?>A' we mus t have e Su Z^(u) — > almost surely by using a familiar Borel- 
Cantelli argument. [The condition pip^ < tjj^ guarantees that the expec- 
tation in the last line above tends to a finite limiting value, hence stays 
bounded over all times s, as can be checked by using formula (17), for ex- 
ample.] 

For all < p < 1 we find pip^ < ippy Considering the graph of we 

quickly see that, for A £ [A(6*),0), taking p as close to 1 as we like gives the 
best rate. For A £ [A m j n , X(9)) we can choose p so that pX = \(6), which gives 
the best rate. 

Recall from Theorem 17 that Z^(oo) > when A £ (A(6>), 0). Then, so far, 
we have proved the following: 

Lemma 22. For every starting law, P x ' y , and for all e > 0, if A £ (A m i n , 0) 
then 
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where 

* . = (c(9), if x mm <x<\(e), 

A ' K, if \(6)<\<0. 

It is clear that this gives the required upper bound of 

iimsup - < A(c A — c A j. 

Now, for any e > 0, if A G (A m i n , A(6>)] then 

e rt e ^Y u (tf+x(x u (t)+c(e)t) > e A(L t +c(e)t)+rf ^ a g 

since we know that Lt := inf{X u (t) : u G iV(f)} satisfies Lt/i — ► —c{0) a.s. 
Otherwise, with AG (A(0),O), 

e rf £ e ^Yu(tr+x(x u (t)+c-t) > e rf Z - (t) ^ oo a . s . 

u£N t 

since here Z A (oo) > a.s. Thus, in all cases, 

u inf !2^£W> A(c ±_ c j), 

t— >oo t 

which completes the proof of Theorem 20. □ 

9. Proof of Theorem 3. Upper bound. The idea for the upper bound 
proof is to overestimate indicator function by exponentials, and then re- 
arrange the expressions to form martingale terms. 

Simply observe that for A G (A m i n ,0), 

iVt(7,[M/t,00))= 1 {X u (t)<-~ / t;Y u (t)> K Vt} 

ueN t 



(64) 



^ n „i>t(Yu(t) -re 2 t)+\(X u (t)+ 1 t) 

uGN t 

< e (S+-re 2 ^++A 7 )t e 4>+Y u (t)2+\X u (t)-E+t 

ueN t 

<e-^ c t~ c x) t Z+(t)e iE ^ + ^- K2 ^)\ 

where Ef = — Ac A . 

Recall from equations (11) and (32) that E^ + A7 — K 2 ip^ has a minimal 
value of A(7, k) achieved when A = A(7,k). Since c{9) is the minimal value 
of c\, Theorem 20 implies that 

(65) limsupt" 1 logZ+(t) < A(4 -c A ) 

t— >oo 
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almost surely for all A G (A m i n ,0). 

In cases where A (7, k) < 0, we can use the optimal value for A, Theorem 20 
and trivially note that ^(7, [n\/t, 00)) is integer valued to deduce that 

E 1 {Y u {t)>K,v r t;Xu(t)<~ 1 t} =0 
u£N t 

eventually, almost surely. Hence, -D(7, k) = —00 almost surely if A(7, k) < 0. 

Otherwise we have A(7, k) > 0, which in fact guarantees that 7 G (0, c{9)\ 
and hence A(7,k) G [A((9),0). Then since 

limsupi -1 logNt(j, [ny/i, 00)) 

< limsupr 1 log(e- A ^- c A >Z+{t)) + (E^ + A7 - k 2 ^ ) 

t— >oo 

we can again make use of Theorem 20 and the minimizing A value, A (7, k), 
to get the bound 

limsupi" 1 log ^(7, [nVt, 00)) < A (7, re) almost surely, 
as desired. 

Notice that, when A(7, re) = 0, the right-hand side of the inequality at 
(64) will tend to infinity (see Corollary 19). Then, on the boundary, we have 
only shown that limsupt -1 log ^(7, [n\/t,oo)) < 0. 

10. Proof of Theorem 1. The spatial growth rate. We first bound the 
spatial growth rate above. Suppose that C C K is Borel-measurable with 
f c e^ y2 <j)(y)dy > 0. Let A G (A min ,0), then 

E hx u{ t)<-^Yumc } < E hYumc } e x{Xu{t)+lt) 

uGN t ueN t 

u£N t 

Recalling equations (8) and (19), we therefore have 

E Hx u (t)<^t;Y u mC} <e A ^Z- (t). 
uGN t 

Now if 7 > c(6>), corresponding to A 7 G (A m i n , A(#)] and having A(7) < 0, 
we know from Theorem 17 that Z^(oo) = almost surely. Then, 

7 > c(0) Hx u {t)<—it;Y a {t)&C} = eventually, a.s. 
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Otherwise, if 7 £ (0, c(6*)), corresponding to A 7 6 (A(6>), 0) and having A(7) > 
0, Theorem 17 tells us that Z^(oo) > almost surely, hence 

limsupt-Mog Hx n (t)<-yt;Y u (t)£C} 

Now we bound the growth rate from below. Let e > be small, A(#) < 

d 2 E~ 

A < 0, and fi = A — e. We recall now that E x is convex so qt£ > and 
7^ > 7a- Then 

E„XX u (t)-E7t^ 
e X H-hx+e)t<X u {t)<-(^-e)t;Y u (t)&C} 

u£N t 

<r c A (-(7A+e)t)--B~t 1l 

^ e 1 {-(7A+e)*<^W<-(7A-£)i;^WeC'} 

u£N t 

— J-*-y\-E~-\e)t fl 
<f C (-^TA-Sr-Ae)t 11 

^ e A 1 {Xu(t)<^(jx~e)t;Y u (t)eC}- 

ueN t 



Then 



rMog £ l { y u ( i)e c } e AX " (i) "^'l { | Xti(i)/t+7A | <£} 

<-A 7A -^-Ae + r 1 log 2 l {Xu(i) <_ (7A _ e)4; y u(i)GC} . 

Letting i — > 00, using Corollary 19 and remembering that for X(9) < A < 
we have Z^(oo) > a.s., we find 

< -A 7A - El - Xe + limmf t~ l log ^ l {x „ (t )<_ (7A _ £)t;Yu(t)eC} 

and as e > can be arbitrarily small we have 

U t ^™ f * _1 log ^ 1 {^W<-7Ai;^(*)eC} > + A 7A- 
«ev t 

Equivalently, 

uminf i _1 log ^ l{x„(t)<- 7 t;y„(t)6C} > ^ 7 + A 77 = A( 7 ) 

ueJVt 

and hence the lim sup and lim inf agree as required. 

We note that these proofs will easily adapt to cover a multi-type branching 
Brownian motion where the types evolve as a finite state Markov chain, 
such as found in [2], where it will also be possible to prove the analogous 
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convergence theorem required when we have a finite type space by adapting 
the proof of Theorem 18 found in [12]. 

In the standard branching Brownian motion case things are even simpler 
to adapt (where, of course, there is no need for any convergence result akin to 
Theorem 18). All the information necessary is contained in the martingales 
J2u£N t ex ^P(^X u (t) — (A 2 /2 + r)t) studied by Neveu [22] and, as first came 
to our attention during discussions with J. Warren, the martingale with 
parameter A can only be capable of "counting" particles near in space 
at large times t, so when this martingale is uniformly integrable particles 
must perpetually be found with the corresponding speed. Of course, in this 
case more precise results, in the spirit of Watanabe [25], also exist. 
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